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In order to observe a signal of possible CP violation in top-quark couplings, we 
have studied energy correlation of the final leptons in e + e~ — > tt — > £ + £~X / £ X 
at future linear colliders. Applying the recently-proposed optimal method, we 
have compared the statistical significances of CP-violation-parameter determi- 
nation using double- and single-lepton distributions. We have found that the 
single-lepton-distribution analysis is more advantageous. 
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The top quark, thanks to its huge mass, is expected to provide us a good 
opportunity to study beyond-the-Standard-Model physics. Indeed, as many au- 
thors pointed [1 — 8] , CP violation in its production and decay could be a useful 
signal for possible non-standard interactions. This is because (i) the CP violation 
in the top-quark couplings induced within the SM is far negligible and (u) a lot 
of information on the top quark is to be transferred to the secondary leptons 
without getting obscured by the hadronization effects. 

In a recent paper, we have investigated CP violation in the tt-pair pro- 
ductions and their subsequent decays at next linear colliders (NLC) ||. We 
have focused there on the single-lepton-energy distributions. In this note, we 
study both the double- and single-lepton-energy distributions in the process 
e + e~ — ► tt — > £ + £~X / £ ± X, and we compare the expected precision of CP- 
violation-parameter determination in each case. For this purpose, we apply the 
recently-proposed optimal procedure ||. 

Let us briefly summarize the main points of this method first. Suppose we 
have a cross section 



where the fi(4>) are known functions of the location in final-state phase space 
and the q are mo del- dependent coefficients. The goal would be to determine 
q's. It can be done by using appropriate weighting functions Wi(4>) such that 
/ Wi(4>)T,((j))d(f) = q. Generally, different choices for Wi(<f>) are possible, but there 
is a unique choice such that the resultant statistical error is minimized. Such 
functions are given by 
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where Xi~ is the inverse matrix of M,-,- which is defined as 




(2) 
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When we take these weighting functions, the statistical uncertainty of q becomes 



Aa = JXu a T /N . 



(3) 



where <jt = J(da/d<j))d<j) and N = L e s<jT is the total number of events, with L e ff 
being the integrated luminosity times efficiency 

In our analyses, we assume that only interactions of the third generation of 
quarks may be affected by beyond-the-Standard-Model physics and that all non- 
standard effects in the production process (e + e~ — > tt) can be represented by the 
photon and Z-boson exchange in the s-channel. The effective jtt and Ztt vertices 
are parameterized in the following form 



v{p t ), 



(4) 



(v — 7 or Z) 

where g is the SU(2) gauge- coupling constant. In principle, there are also four- 
Fermi operators which may contribute to the process of tt production. However, 



as it has been verified in Ref . |L(J , their net effect is equivalent to a modification 
of A v and B v . Therefore, without loosing generality we may restrict ourself to 
the vertex corrections only. 

For the on-shell W, we will adopt the following parameterization of the tbW 
vertex: 



where P, 
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u(Pt), (5) 
v(pb), (6) 
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'l=F7 5 )/2, V t b is the (tb) element of the Kobayashi-Maskawa matrix 



and k is W*s momentum. 

Using the above parameterization, applying the narrow-width approximation 
for the decaying intermediate particles, and assuming that the Standard-Model 
contribution dominates the CP-conserving part, we get the following normalized 
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double- and single-lepton-energy distributions of the reduced lepton energy V= 



2EJ (1 — /?)/(! + P)/m t , E being the energy of £^ in the e + e cm. system, and 



(3 = ^l-Am^/s: 
Double distribution 



1 d 2 a 3 

22cifi(x,x), (7) 



a dx dx , | 



where x and x are for £ + and £ respectively, 



and 



Cl = l, c 2 =£, c 3 = ^Retf* - jj 



/i(ar, x) = f(x)f(x) + rf g(x)g(x) + r] [ f(x)g(x) + g(x)f(x) ], 
f 2 (x, x) = f(x)g(x) - g(x)f(x), 

f 3 (x, x) = 5f(x)f(x) - f(x)6f(x) + rf [ 5g(x)g(x) - g(x)6g(x) ] 
+rj [ 5f(x)g(x) - f(x)5g(x) + 5g(x)f(x) - g(x)5f(x) }. 



Single Distribution 



1 <] " E*), (8) 



a ± dx , i 
where ± corresponds to £ , 

cf = l, cf = Tt, 4 = Mf?), c 3 - = Re(/ 2 L ) 

and 

h(x) = f(x) +r] g(x), f 2 ( x )=g(x), f 3 (x) = 5f(x) + rj Sg(x). 

Since all the functions and parameters in these formulas are to be found in Refs. |7|, 
IH, we only remind here the normalization of f(x), 5f(x), g(x) and Sg(x): 

f(x)dx = 1, / 5f(x)dx = / g(x)dx = / Sg(x)dx = 0. (9) 
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7], v( and £ are numerically given at y/s = 500 GeV as 

T) = 0.2021, if = 1.3034, £ = -1.0572 Re(£> 7 ) - 0.1771 Re(D z ) 

for the SM parameters sin 2 6 W = 0.2325, M w = 80.26 GeV, M z = 91.1884 GeV, 
P z = 2.4963 GeV and m t = 180 GeV. 

In Eqs.(0,D, CP is violated by non- vanishing £ and/or Re(f^ — / 2 ) terms.0 
First, let us discuss how to observe a combined signal of CP violation emerging 
via both of these parameters. The energy-spectrum asymmetry a(x) defined as 

_ da~ /dx — da + / dx 
da~ I dx + da + / dx 

has been found as a useful measure of CP violation via £ @, 0- In Ref. || we have 
computed a(x) for the case where both £ and Re(/^ — f£) terms exist. Practically 
however, measuring differential asymmetries like a(x) is a challenging task since 
they are not integrated and therefore expected statistics cannot be high. For this 
reason, we shall discuss another observable here. 
A possible asymmetry would be for instance 

dxdx- — : / / dxdx- 



\ J Jx<x uxax J Jx>x 

A u = —— jti — -w— ■ l 1U J 



•. r dxdx J Jx>s dxdx 

d^cr r r d^cx 
dxdx- — — + / / dxdx- — — 
x<x dxdx J Jx>x dxdx 

For our SM parameters, it becomes 

A u = 0.3638 Re(£> 7 ) + 0.0609 Re(D z ) + 0.3089 Re(/ 2 fi - f^) 

= -0.3441 £ + 0.3089 Re(/ 2 * - f 2 L ). (11) 

For Re(D 7 ) = Re(D z ) = Re(/ 2 i? ) = -Re(/ 2 L ) = 0.2, e.g., we have 

A ff = 0.2085 



^In the present note, t, t and W are assumed to be on their mass shell since we are adopting 
the narrow-width approximation for them, and the contribution from the imaginary part of the 
Z propagator is also negligible since s is much larger than M§. Therefore we do not have to 
consider CP- violating effects triggered by the interference of the propagators of those unstable 
particles with any other non-standard terms Jill] . 
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and its statistical error is estimated to be 

AA U = - A 2 U )/N U = 0.9780/ 

Since a eS ^ t t — 0.60 pb for y/s = 500 GeV, the expected number of events is 
Nu = 600 euLBj, where e« stands for the £ + £~ tagging efficiency (= ej ; q is 
the single-lepton-detection efficiency), L is the integrated luminosity in fb _1 unit, 
and Be(~ 0.22) is the leptonic branching ratio for t. Consequently we obtain the 
following result for the error 

AAu = 0.1815/ V /^L, (12) 

and thereby we are able to compute the statistical significance of the asymmetry 
observation Nsd = \Age\/AAee. 

In Fig.l we present lines of constant Nsd as functions of Re(_D 7 ) = Ke(Dz) 
and Re(/^ — f£) for L = 50 fb _1 and e« = 0.5 (which mean Nu = 726). Two 
solid lines, dashed lines and dotted lines are determined by 

| 0.4247 Re(ZV) + °- 3089 Mf? " J2) I = N SD /^N 2 SD + N U 

for Nsd =1) 2 and 3 respectively. We can confirm An to be non-zero at la, 2a 
and 3a level when the parameters are outside the corresponding lines. It can be 
seen that we have good chances for observing the effect at future NLC unless 
there is a conspiracy cancellation between those parameters. Table |IJ shows the 
a/s dependence of Nsd for the same e«L. 

In order to discover the mechanism of CP violation, however, it is indispens- 
able to separate the parameter in the top-quark production and that in the 
decay (Re(/^ — /2O). We shall apply the optimal procedure of Ref.|| to the dou- 
ble distribution first. Using the functions in Eq.(^), we may calculate elements 
of the matrix M and X defined in Eqs.(|I], |2[): 

M11 = 1, M 12 = M X3 = 0, M 22 = 0.2070, M 23 = -0.3368, M 33 = 0.6049 

* 2 We use £ instead of Re(D 7 z) as a basic parameter when we discuss parameter measure- 
ments, since £ is directly related to the distributions Eqs.(^J) and (||). 
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Figure 1: We can confirm the asymmetry Au to be non-zero at la, la and 3a 
level when the parameters Re(-D 7 ^) and Re(/^ — f 2 ) are outside the two solid 
lines, dashed lines and dotted lines respectively. 

and 

X n = 1, Xn = X 13 = 0, X 22 = 51.3389, X 23 = 28.5825, X 33 = 17.5662. 
This means the parameters are measured with errors of 

A£ = 7.1651/ yfifo, AReiff - f})(= 2^X 22 /N U ) = 8.3824/ fi^. (13) 

Next we shall consider what we can gain from the single distribution. We 
» 3 Note that a T in Eq.(|) is unity in our case since we are using normalized distributions. 
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y/s (GeV) 


500 


600 


700 


800 


900 


1000 


0~ee^tt (P b ) 


0.60 


0.44 


0.33 


0.25 


0.20 


0.16 


P = 0.1 


2.8 


2.5 


2.3 


2.0 


1.8 


1.7 




(0.1043) 


(0.1097) 


(0.1132) 


(0.1155) 


(0.1171) 


(0.1183) 


P = 0.2 


5.7 


5.2 


4.6 


4.1 


3.7 


3.4 




(0.2085) 


(0.2195) 


(0.2263) 


(0.2309) 


(0.2342) 


(0.2365) 


P = 0.3 


8.9 


8.0 


7.2 


6.4 


5.8 


5.3 




(0.3127) 


(0.3292) 


(0.3395) 


(0.3464) 


(0.3513) 


(0.3548) 


P = 0.4 


12.4 


11.3 


10.1 


9.1 


8.2 


7.5 




(0.4170) 


(0.4389) 


(0.4527) 


(0.4619) 


(0.4683) 


(0.4730) 



Table 1: Energy dependence of the statistical significance N$d of A(± measure- 
ment for CP-violating parameters Re(P> 7 ) = Re(D z ) = Re(f^) = — Re(/^)(= 
P) =0.1, 0.2, 0.3 and 0.4. The numbers below N$d (those in the parentheses) are 
for the asymmetry 



have from Eq.(§) 

Mh = 1, M 12 = Mi 3 = 0, M 22 = 0.08 98, M 23 = 0.1499, M 33 = 0.2699 
and 

1, X 12 = X 13 = 0, X 22 = 151.9915, X 23 = -84.4279, X 33 = 50.6035. 



X 



li 



Therefore we get A£ = 12.3285/ y/Ni and ARe(f?) = 7.1136/ v / iV7 from the £ + 
distribution, and analogous for At; and Z\Re(/ 2 L ) from the l~ distribution. Since 
these two distributions are statistically independent, we can combine them as 



Ai = 8.7176/ 'JN E 



ARe(f 2 R 



fl) 



10.0601/ J N e . 



(14) 



It is premature to conclude from Eqs. ([l3D and (14) that we get a better 
precision in the analysis with the double distribution. As it could be observed in 
the numerators in Eqs. (|T3"|, [14]), we lose some information when integrating the 
double distribution on one variable. However, the size of the expected uncertainty 
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depends also on the number of events. That is, Nu is suppressed by the extra 
factor comparing to Ng. This suppression is crucial even if we could achieve 
tp — \. For N pairs of tt and q = 1 we obtain 

Ai = 32.5686/^, ARe(f* - ft) = 38.1018/ViV 

from the double distribution, while 

A£ = 18.5859/^, ARe(f? - f%) = 21.4484/ViV 

from the single distribution.^ Therefore we may say that the single-lepton-distri- 
bution analysis is more advantageous for measuring the parameters individually. 

In summary, we have studied how to observe possible CP violation in e + e~ — > 
tt — > £ + £~X and £^X at NLC For this purpose, CP-violating distributions of the 
fmal-lepton energies are very useful. Using these quantities, we introduced a new 
asymmetry Au in Eq.(|II]), which was shown to be effective. Then, applying the 
optimal procedure 0, we computed the statistical significances of CP- violation- 
parameter determination in analyses with the double- and single-lepton-energy 
distributions. Taking into account the size of the leptonic branching ratio of the 
top quark and its detection efficiency, we conclude that the use of the single-lepton 
distribution is more advantageous to determine each CP-violation parameter sep- 
arately. 
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